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It is shown how continuous products of random transformations con- 
strained by a generic group structure can be studied by using Iwasawa's 
decomposition into "angular", "diagonal" and "shear" degrees of freedom. In 
the case of a Gaussian process a set of variables, adapted to the Iwasawa 
decomposition and still having a Gaussian distribution, is introduced and 
used to compute the statistics of the finite-time Lyapunov spectrum of the 
process. The variables also allow to show the exponential freezing of the 
"shear" degrees of freedom, which contain information about the Lyapunov 
eigenvectors. 
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1. Introduction 

Lyapunov exponents represent the natural way to measure the rapidity of divergence of 
a set of infinitesimally close points carried along by a chaotic or random flow [31 1141 IT]. 
The simplest example is given when a smooth TV-dimensional dynamical system of the 
form 

x = F (x) 

is linearized around a trajectory x (t) giving the linear equation y = X (x) y for the 
separation y (t) between two infinitesimally close orbits. The generator of the tan- 
gent dynamics, Xij — d{F/dxj, is known as the stability matrix of the system. The 
evolution of the initial separation y (0) is formally given by the expression 

y(t)=ff(t)y(Q) 
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where g (t) is the time-ordered product 



g (t) = T exp 



X (t') At' 



The N time-dependent Lyapunov exponents Ai (t) > ■ ■ ■ > 
eigenvalues of the symmetric matrix (2t) _1 log [g^ (t) g (t)] 



(1) 

Xn (t) are defined as the 
and measure the rate of 



N). Oseledec (S3 proved 
oo limits of the (<) and of 



growth of infinitesimal fc-dimensional volumes (fc — 1 
under very general conditions the existence of the t - 
the eigenvectors it (t) of g' (t) g (t). 

A common ansatz for strongly chaotic systems assumes that the stability matrices 
X (t) decorrelate in a characteristic time, shorter with respect to other time scales 
of the system, and can therefore be treated as a succession of independent linear 
transformations ^3j. Then the separation y (t) evolves according to the stochastic 
equation 



y = *(*)y (2) 

After carefully specifying the statistics of the random process X (t) , and choosing in 
particular one of the standard (Ito or Stratonovich) regularization for (2J , the averages 
of functionals / [y (t)] over the random realizations y (t) can be written in the form of 
path integrals |34j : 

(/[y (*)]> = JvXf{y(t)]e- s W (3 ) 

The problem of computing averages of the form © is quite general, and appears both 
in the context of quantum and classical systems. 

In the theory of quasi one-dimensional disordered wires P] the study of low-temperature 
conductance fluctuations in small metallic samples is reduced to the computation of 
the statistics of g^ (L)g(L), where g (L) is the 2n x 2n transfer matrix for a wire 
of length L and n is the number of transverse modes at the Fermi level pi]- The 
transfer matrix g (L) can be naturally written in the form as an ordered product 
of transfer matrices for infinitesimal portions of wire. Here the role of time is played 
by the length L of the sample, and the role of the Lyapunov exponents by (inverse) 
localization lengths. Ensemble averaging implies the computation of path integrals of 
the form , where the infinitesimal transfer matrices X are constrained by symmetry 
considerations. In the presence of time reversal symmetry the X matrices belong to 
the sp (n, R) algebra, if time reversal symmetry is broken they belong to su (n, n), and 
in the presence of spin-orbit scattering to so* (4n) |25j . 

In the theory of passive scalar transport in n-dimensional space the computation 
of many statistical quantities, such as the exact p.d.f.'s of the scalar ^l] EE GDI E] and 
of its gradients ^OlE], can be reduced to the form (JHJ, where X is a matrix of velocity 
gradients of the carrying fluid. The matrices X are the infinitesimal generators of the 
dynamics of an incompressible fluid element and belong to the si (n, R) algebra. 
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In the study of deterministic chaos, products of random matrices have been used to 
mimic the chaoticity of dynamical systems [Oj. I n the case of n-dimensional Hamil- 
tonian systems one has products of random symplectic maps [HJ EB] . This formally 
leads to expressions of the form <EJ where the matrices X belong to the sp (n, R) al- 
gebra. Similar models have been recently used to study the hydrodynamic modes in 
the Lyapunov spectrum of extended dynamical systems 

In many of these cases analytic results were obtained. Such results were made 
possible by the presence of three main ingredients: i) Markovianity 1 , ii) Gaussianity, 
and Hi) the presence of a group structure in the set of transformations and of a 
statistics compatible with it (i.e., a statistics that can be expressed in terms of the 
intrinsic geometry of the transformation group) . 

While the physical support for the first two assumptions is in some cases question- 
able, the existence of a group structure lays on stronger foundations, since it directly 
reflects the symmetries of the underlying physical problem. 

In the case of quasi one-dimensional wires, results were obtained studying the 
Fokker-Planck equation for the transmission eigenvalues as a diffusion equation over 
Riemannian manifolds ^JEj- A similar approach has been also applied to the case of 
passive scalar transport [ZJ. On the other hand, functional integral representations of 
the form <P} for the probability densities can be exploited directly j.'ill |H1 111! 1131 1101 
1211 112j and can provide an intuitive description of the statistics of the Lagrangian tra- 
jectories described by Eq. J2j. Path integral representations also provide the ground 
for several approximations, like instanton calculus and perturbative expansions. Ac- 
tually, the two approaches are complementary and one should switch from one to the 
other depending on convenience. However, the Lagrangian theory has been worked out 
mainly on specific examples, often relying on particular features of the given algebras 
|20l I21j . and not covering the most general case. The purpose of this work is to present 
this theory in a sufficiently general setting. It turns out that a formal exposition not 
only allows to cover many different cases in a unified way, but also unveils the simple 
geometric meaning of some formula? and presents a very simple and appealing scenario, 
which can hopefully serve as a basis for further investigations. 

The starting point of the Lagrangian theory is a convenient decomposition of the 
time-ordered exponential (QJ. Let us remind here some facts about the standard 
method jH] for the numerical computations of Lyapunov exponents. An orthonormal 
set of k vectors is evolved for a time t and a Gram-Schmidt orthogonalization pro- 
cedure is applied to the evolved vectors in order to obtain a new orthonormal set. 
The procedure is iterated, and by performing an average of the logarithmic fc-volume 
increase over many steps one obtains the approximated sum of the first k Lyapunov 
exponents. When the evolution operator g (t) belongs to the SL (n, R) group it is easy 
to see j2H] that the Gram-Schmidt orthogonalization procedure is equivalent to the 
decomposition of the evolution operator in the product or three consecutive transfor- 
mations, corresponding to a shear, a dilation, and a rotation: the rotation carries the 



Although results were obtained also with time-correlated statistics, see e.g. Refs. [HlllS] . 
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old orthonormal set in the new one, the dilation is responsible for the volume increase, 
and the shear part inherits from the Gram-Schmidt procedure a triangular form. From 
an algebraic point of view, the group element g (t) is decomposed into the product of 
a compact, Abelian and nilpotent part. Iwasawa [26_ showed that a similar decompo- 
sition can be defined for a large class of groups and is a smooth mapping (see Ref. jJH] 
and App. EJ. 

The fact that emerged from a series of works (221 EH E0 1201 IS] is that starting 
from normally distributed, centered random variables X it is possible in certain cases 
to introduce new variables, adapted to the Iwasawa decomposition 2 , which are still 
normally distributed. Here it is shown that the new variables can be introduced in 
the quite general case where the random variables X belong to the algebra of a real 
semisimple Lie group, and their covariance has a form compatible with the group 
structure (c/. Eq. ©). It turns out that the definition of the new variables (the 
{A, M} variables defined by Eqs. Il.'>l2.ili ) is most easily given in terms of the group 
structure of the set of transformations, and in particular of the "adjoint operator" 
M(g)X = gXg-\ 

The fact that the new variables are normally distributed is non trivial, since the 
Iwasawa transformation is nonlinear, and it is of course quite relevant in practical 
calculations. In particular, it allows the rapid computation of the Lyapunov spectrum, 
which is encoded in the statistics of the Abelian part of the Iwasawa decomposition. 
Although the original variables X are centered, the Abelian part acquires a positive 
mean value thanks to the non-trivial Jacobian factor associated with the variable 
transformation (c/. (l35l) L This mean value is precisely the finite-time Lyapunov 
spectrum, which can be expressed explicitly in terms of the structure constants of the 
algebra of the X variables (cf. Eqs. QHESl)- 

The new variables however contain more information than just the statistics of 
Lyapunov exponents. For instance, Eqs. H42I38I37I) prove in a quite general setting 
the freezing of the "shear" (nilpotent) degrees of freedom which was used in Refs. 
[lOl 121] to compute the statistics of passive scalar dissipation. The nilpotent degrees 
of freedom actually encode information about the Lyapunov eigenvectors f& (t). For 
instance, in the SL (n, R) case it is easy to see that the ffc (t) can be obtained (up 
to exponentially small terms) by Gram-Schmidt orthogonalization of the row vectors 
of the (nilpotent) shear matrix. The exponential freezing of the shear part therefore 
implies the exponentially rapid convergence of the eigenvectors f& (t) to their limits for 
t — > 00, thus providing a particularly clear realization of Oseledec theorem. 

The work is organized as follows. In Sec. I define the random process X (t) by 
assigning its covariance in terms of intrinsic geometric features of the corresponding 
algebra. This is equivalent to assigning a Feynman-Kac measure on the space of paths 
t — » X (t). The multiplicative process g (t) is then defined as a time-ordered expo- 
nential through the appropriate (Stratonovich) regularization. In Sees. I3I4I5I the new 



2 Or to the Gauss decomposition, which was considered in Refs. [291 il3j respectively for the su (2) 
and su (n) cases. 
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variables are defined through the following steps: i) performing the Iwasawa decom- 
position of the process g(t); ii) introducing intermediate variables, adapted to the 
Iwasawa decomposition; Hi) rewriting the Gibbs weight exp (— S) in terms of the new 
variables; iv) computing the Jacobian factor of the variable transformation, i.e., the 
transformation rule for the "volume form" VX; v) performing a Lie algebra transfor- 
mation that leads to new, normally distributed variables. Although the intermediate 
steps are technical, the final results are simple and are summarized by formula? (|36l 
14011 . In Sec. El the new variables are used to compute the statistics of the finite time 
Lyapunov spectrum of a generic linear representation of the process g (t). In Sec. 0the 
cases of the si (n, R) and sp (n, R) algebras are worked out in detail. A brief summary 
of results of Lie algebra theory are given in the Appendix. 

2. Random walks on Lie groups 

The computational scheme outlined in the Introduction can be implemented in the 
quite general case where the set of transformations g (t) form a real semisimple Lie 
group G (this comprises the case of complex groups, as soon as they are seen as real 
groups endowed with a complex structure) . The infinitesimal generators X (t) then 
lie in the corresponding algebra g. Each real semisimple Lie algebra is characterized 
by two natural structures: a canonical quadratic form B {X, Y) known as the Killing 
form, and an involution 9 such that B {X, 9Y) becomes strictly negative definite (see 
Appendix El . For instance, in the case g = si (n, R) (the algebra of real n x n real 
matrices with null trace) one has B (X, Y) = 2nTrXY, 9X = — X tT and the effect of 
9 is that of changing the sign of the symmetric part of X. 

In order to define a Gaussian probability density on g one needs a positive definite 
quadratic form. I will consider here the general linear combination 

F (X, Y) = fiB (X, Y) + vB (X, 9Y) (4) 

where /i, v are real parameters such that F be positive definite (see itlfill below). 
Let X (t) be the Gaussian process on [0, T] with values in g, zero mean and covariance 

{X(s)®X(t)) = 5{t-s)-C (5) 
where C = F^ 1 , i.e. the random process defined by the Feynman-Kac measure 

T 

l[dX(t) (6) 

t=o 

where dX is a Euclidean measure on g and K is an (infinite) normalization constant. 
The continuous process X (t) can be seen jlfJI I22j as the limit of the discrete process 



dfi (X) = K exp 



F(X (t),X(t))dt 
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Xj, where Xj = X (tj), tj = je, e = T/N, defined by the probability measure 



d/ijv (X) = K n exp 



1 N 

--^zFiX^Xj) 



3 = 1 



N 



(7) 



as N — > oo . 

Let g (t) be the time-ordered exponential defined by 



g (t) = T exp X (t) dr) = ^lim 



fjN 



where gj is defined recursively by 

9] = exp (eX^gj^x, 



go = e (the identity) 



(8) 



This discretization prescription is compatible with the group structure, and is of 
Stratonovich type. As a matter of fact, if G is a linear group (and, in general, on 
the universal enveloping algebra of g) one can compute 



1 

2~c 
1 

1 

Ye 

X s + O (e 2 X3) 



(9 j -9j-i) {gj 1 + .97-1 ) 

(gjgj\ - gj-igj 1 ) 

[exp (eXj) - exp (-eXj)] 



obtaining 



9j~9j-i 9j ■ .'/., : 



where O (e 2 ^ 3 ) denotes a series of terms having order at least 2 in e and at least 3 
in Xj. Such terms can be neglected in averages when N — > oo, as usual in the theory 
of Feynman-Kac integrals pH], since (Xj <g) Xi) = e~ 1 Sji ■ C ~ 1/e. In other words, in 
the limit N — > oo one can consider 

^(t) = .g(t)g- 1 (<) 

understanding that the midpoint regularization JHJ has to be used throughout. 



3. Decomposition of the trajectories 

Following the general scheme outlined in the Introduction, let us perform the Iwasawa 
decomposition 115615711 of the group element gj corresponding to the discrete time j in 
the product of a compact, Abelian and nilpotent part: 

9j = (10) 
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where gj is defined recursively as in (jHJ and kj , aj , rij are functions of gj obtained 
by inverting the diffeomorphism ll58l) . Moreover, let us define Kj, Aj, Nj through the 
relations 

kj = cxp (eKj) kj-i, aj = cxp {(-Aj) o,-_i, rij = cxp (eNj) (11) 

The variables K, A, N are an intermediate set of new variables, adapted to the Iwasawa 
decomposition 111 Oil . In order to obtain the explicit form of the variable transformation 
X — » {K,A,N}, one substitutes {Ell m ® and uses l(TT)l : 

exp(eXj) = gjgj\ 

= exp (eKj) ■ Ad (kj-i) exp (eAj) ■ Ad (fcj_iaj_i) exp (eAj) 

where Ad (<?) A = gXg~ x . Using © one obtains: 

Xj = l[exp( e A J )-exp(-eA J )] + 0( e 2 A J 3 ) 

= K 3 + Ad (fcj_i) Aj + Ad (fej_ 1 a i _ 1 ) A,- (12) 
+ | [Aj, Ad (fej_x) A,-] + | [Aj, Ad (fcj-iaj-i) A,-] 

+ | [Ad (kj-i)Aj,Ad (fcj-iOj-i) Aj] + O (e 2 A|) 

where O (e 2 A|) denotes a series of terms of at least second order in e and at least 
third order in Xj, Kj, Aj, Nj. Using 

Ad (exp (^x))Y = Y+ e -[X,Y]+0{e 2 X 2 Y) 

this can be rewritten as 

Aj = Kj + Ad (kj-A A 3 + Ad {kj-xCLj-^ N 3 + O (e 2 X]) (13) 

with kj-i = exp (§Kj) fej— i, Sj-i = exp (f-Aj) Oj-i. 

Formula l(T3ll extends the Iwasawa decomposition on the Lie algebra g to a natural 
decomposition 

(K j ,A j ,N j )^X j 

of the random process Xj , which corresponds to the decomposition ifTfljl at the group 
level. Observe that 

= %-i {KjiKj-Xi ■ ■ ■ > Ki) , a-j-i = C4j>4j-i> . . . ,A X ) 

and that the appearance of kj_i and means that at the continuum level the 

midpoint regularization 10 applies also to the group terms arising from the Iwasawa 
decomposition. 
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4. Decomposition of the Gibbs weight 



Relation ltl3ll can now be used to express the Gibbs weight defined through J3J in 
terms of the new variables. Using the orthogonality with respect to B of a to both k 
and n, of n with itself, and properties ^531541551 , a simple computation shows that 



with 



FiX^Xi) = { l i + v)B{K j -K° j ,K j -K^) 
+ (fi-v)B(A j ,A j ) 

-^-HB{kd{aU)N^e N] ) 
Ad (kj-iaj-ij N 3 +0Ad (kj^aj-ij Nj 



(14) 



K° = -- 

J 2 



Recall ing tha t B(X,X) > on p, B{X,X) < on k and B (X, OX) < on g, and 
using H54I55H . it is clear that F (X, X) > if 



)jl > v and fi < —v 



(15) 



The last two rows of 114t can be written more explicitly by expanding over the basis 
Y a defined in (ESI and a basis Hi of a: 



aeP 



Then using 



with 



one gets 



ad (Aj)Y a = a(Aj)Y a 
Ad (<ij-_ i) Y a = exp [a (exp -1 Y a 



exp 1 a^i =J2 eA i+ ~ A j 



i=i 



B(A J: A J ) = 
Y,A k j A]B{H k ,H l ) 



k.l 



B (Ad (a^N^eNj) = 

N?N? exp [2a (exp -1 a^)] B (Y a ,6Y p ) 

a,peP 



(16) 



(17) 
(18) 



(19) 



(20) 



(21) 
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Here the matrix B (Hk, Hi) is directly related to the Cartan matrix which characterizes 
the Lie algebra g. A list of Cartan matrices for all semisimple Lie algebras can be 
found in Ref. |23]- The matrix B(Y a ,6Yp) is non-degenerate, it is diagonal if the 
multiplicities l|59fl are all 1 and is related to the way the Cartan involution 9 permutes 
the roots otherwise |23] . 

It is often necessary to compute averages of functionals F ({aj}) which do not de- 
pend on kj, rij. In this case, the form of l(T4T) shows that integration over YljLi dKj 

produces only a constant factor. Integration over Y[f =1 ^Y? produces instead a factor 
proportional to 



D 



-1/2 : _ 



N 



N 

n n exp \-~ a ( ex p _i 

j=l a£P 
N 

I [ exp [-2p (exp -1 a,-_i)] 
3=1 

n r j-i 

-2^ep(^)-ep(^) 



(22) 



J|exp 

3=1 



where notation 116511 was used. 

Expression 11411 takes a simpler form if one introduces the "rotated" variables 



M 



(23) 



The variable change N - 
d({Mj}) 



M introduces however a Jacobian factor 

N 



IJdettAdfa-OlJ 

3=1 
N 

J^J exp [Tr ad (exp -1 2j- x) | n 

3-1 

2j2^P(Ai)+ep(A 



3 = 1 
N 

3 = 1 
N 

n° x p 

3 = 1 



(exp 1 a j -i) 



1=1 



D 



1/2 
N 



(24) 



Integration over the variables Mj would than produce a constant (independent from 
Aj) factor. Moreover, factor 12411 will be seen to cancel exactly a corresponding factor 
emerging in the global Jacobian j30jl . 
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5. Decomposition of the volume form 



Along with the variables Xj and Kj , Aj , Nj , which belong to the algebra g of infinites- 
imal transformations, one can consider the "integrated" variables gj and kj,aj,nj de- 
fined in JHJ and lllOt . All together, one has the following commutative diagram of 
variable transformations: 

{.'/.-I — ► I 1 ' , ■",■".;} 

I i (25) 

{Xj} — {Kj,Aj,Nj} 

where the variables of the upper row belong to the group G and those of the lower row 
to the corresponding algebra g. The relation between the two set of variables can be 
found by varying lfl2l : 



SXj = SKj + Ad (kj-i) 5A 3 + Ad (kj -iaj-i) SNj 

+-{ [SKj, Ad (kj-\)Aj\ + [K 3 ,Ad (kj-i)5Aj] 

+ [SKj, Ad fa-m-i) Nj] + [Kj, Ad (kj-iaj-i) SNj] (26) 
+ [Ad (kj-i) SAj, Ad (kj-xdj-i) Nj] 

+ [Ad (kj-i)Aj,Ad (kj-xdj-i) SNj] | 

+0 (e 2 X]5Xj) + f (SXj^SXj^, SXx) 

The variation SXj € Tg can be identified with the 1-form SXj ( • ) = F (SXj , ■). In 
(12611 . O (e 2 XjSXj) denotes a 1-form in SKj, SAj, SNj whose coefficients are at least 
of second order in e and in Kj, Aj, Nj, while / (SXj-i, SXj_2, ■ ■ ■ , SXi) denotes a 
1-form in 5K h SA h SN t for I < j. 

On the other hand, the measure dfi (X) can be re-expressed in terms of the Haar 
measure on G, by computing the variation 



S 9j9j 



eSXj + -(SX j X j + X j SX j ) 
-f{SX j - 1 ,...,SX 1 ) 



eSXi 



(SXj X, + Xj SXj 



-0(e 3 X]SXj)+f(SXj^, 



eSX j + -[X j ,5Xj} + 0(e*XjSXj 



O (e 3 X 2 SXj) 



(1 - *Xj) 
SXi) 

-f(5Xj_ x , 



■9j-i9j 



,5Xx) 
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so that 



Ug 3 gj x = Ad(exp(|x,)) SXj 



+0 (e 2 X]8X J ) + f (SX^ U 5X{) 

In the continuum limit, terms O (e 2 XjSXj) contribute to the measure with a constant 
factor which ensures the correct normalization 3 . One then has 

N dimg N dim G 

A A A *w97 l 
j=i 1=1 j=\ 1=1 

where the 8\X span T g, Sig g _1 span T e G ~ T g, and one uses the facts that 

1. the / (SXj_i, . . . , SXi) terms cancel in the exterior product with SXj terms; 

2. semisimple Lie groups are unimodular, i.e., |det Ad(<?)| = 1 for every g € G. 

So: 

N N 

ll' V * 1I<1</ 

j=l J=l 

where dg denotes the Haar measure on G. In other words, the variable transformation 
{Xj} — > {gj} has trivial Jacobian and d/i (X) = d/i (g), where 



d/j, (g) = K' exp 



~J T F{gg-\gg- l )dt 



II d 5W 



To the two horizontal arrows of diagram (l25l) there correspond the Jacobians defined 



by 



N 



II d 9o = J N({ k j, a j , n 3 }) dkj da 3 drij 

3=1 3=1 
W JV 

II IV = J%({K J .A 3 ,N J })l[dK J dA J dN J 

j=i i=i 



(27) 



(28) 



The variable transformations corresponding to the two vertical arrows have trivial 
Jacobians, so 

,a 3 ,n 3 }) (29) 
3 One can treat terms ~ t 2 X^dXj as a perturbation in integrals over the measure J7J and use Wick's 



theorem to compute their contribution. They give a contribution ( e X. 



e, but their product 



with any power of eXj higher than is negligible in the continuum limit. The contributions ~ e 
sum up to give a constant factor. 
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The product of Haar measures on both sides of ffijl is invariant under the global 
transformations 

9j -» k R gj, gj -> g 3 n L , with k n e K, n h e N 
which correspond to 

kj — * fcR%, rij — > n^riL 

By induction over TV, this shows that does not actually depend on {kj,nj}. From 
il29)l then follows that J% does not depend on {Kj,Nj}. 

The actual value of the Jacobian can now be computed either at the group or at the 
algebra level. At the group level, using Il(i4t . one gets 

TV 

= n cxp [ 2 /°( ex p~ ia j)] 

3=1 



N 

n° x p 

3=1 


. 1=1 








N 

n° x p 

3=1 


' 3-1 

. i=l 


+ ep (vlj) 


• exp 


3=1 



= < 2 -J^' loc (30) 

The decomposition of into an ultralocal part D 1 ^ 2 and a local part J^' loc is such 
that the ultralocal part cancels exactly the factor l|22)l which appears when integrating 
functionals of the form F ({a? }) over 

Ilf=i dK rUf=i dN i- The ultralocal factor D X J 2 
is also exactly canceled when passing to the "rotated" variables Mj defined in Eq. H23I1 . 

The result Q30I1 can be checked by repeating the computation at the algebra level. 
Recalling that Jfj cannot depend on Kj, Nj, we can set Kj = 0, Nj = in <l26l) . 
getting 



SXj = 5Kj + SAj + Ad (aj-i) SNj 

+Z[8K j ,A J ] + Z[A j ,M(a j -i)6N j ] 

+0 (e 2 X]5Xj) + f {SXj_ u SXi) (31) 
= Ad (exp (-|^)) SKj + SA 3 + Ad (exp (|a,-) 0i _i) ^ 
+0 (JxjSXj) + f {SXj_ u 5Xx) 
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Using ijfiTJl and llfi2l) to expand Nj = YlaeP ^NfYa, one gets 

Ad (exp 0i _i) 6Nj 



aeP 



J'-l 



(=1 



(32) 



Using (HU and © to expand SKj = 8K° + ]T QeP 8Kf (Y a + 0Y a ), where 5K$ E m, 
and recalling (Ififil) . one gets 



ad (A^&K^ 



^a(A,)<5J^(y Q -0y Q ) 

5] a (A,) SKJ [2Y a - (Y a + 6Y a )] 

aeP 

N° -^2a(A j )SKJ(Y a + eY a ) 



a£P 



with N$ e n, and 



(l - |ad (A,)) ^ = Yl ex P [ a SK i ( y « + 



a6P 

_i N o + ^A]S Kj ) 



(33) 



Now using (I31I32I33B and the decomposition ll56l) . observing that the terms / (6Xj-i, . . . , SXi) 
cancel in the exterior product with SXj (the transformation {Xj} — > {i^j, A,-, iVj} is 
triangular in the index j), and that O [e 2 XjSXjj terms contribute only to a normal- 
ization factor (see footnote 01, one gets 



JY 



j n ^ n n cxp 

j=l aeP 



O'-l 



a 



J2 tA i + n A i 



= n exp 

3=1 



3-1 



2^ep(A0 + ep(^) 



z=i 

,loc 



exp 



■ exp 



N 

i=i 



(34) 



as expected. 

Taking into account 1-iOII , H34I1 and i!24|) it is clear that the Jacobian of the transfor- 
mation 

is proportional to the simple local factor 



tK Joe 

— exp 



N 



(35) 
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6. Statistics of Lyapunov exponents 

Classical groups are groups of transformations of R™ or C™ . Let ( , ) be a scalar 
product on R™ or a hermitian product on C™. Let v be a vector of R™ or C™, gw the 
action of the group element g on v. Let g (t) be a random walk on a (real or complex) 
classical group G. Then 

|| 3 (t)v|| 2 = (. 9 (t)v, 5 (t)v) = (.gt (t)fl(t)v,v) 

so that the rate of growth of \\g (t) v|| is related to the highest eigenvalue of g' (t) g (t). 
More generally, the rate of growth of the parallelepiped generated by I vectors vi, . . . ,vj, 
is related to the sum of the first I eigenvalues exp (2Ai (t) t) , . . . , exp (2A; (t) t) of the 
hermitian matrix g' (t) g (t) |14| : 

Ax (t) + • • • + A, (t) = i log VolCgW viA-Agfflv,) 
1W w t B Vol(vi A-"Avi) 

The eigenvalues Ax (t) , . . . , Xi (t) are the (time dependent) Lyapunov exponents of the 
process g (t). 

For the classical algebras the Cartan involution 9 is just 6Z = —Z\ which expo- 
nentiates to the group automorphism 6 defined by 65 = (gi) . Let T be a finite 
dimensional representation of G such that (fifij) holds. In general, the (time dependent) 
Lyapunov exponents of the process g (i) can be seen as the eigenvalues of 



(2ty 1 log [THg)T(g)] = - (2ty 1 logT (g^Qg) = {2t)~ l log T (Gr^ 1 • a 
or equivalently as the eigenvalues of 

(2t) _1 log [T (a 2 ) T (nOn" 1 )] 



2 n) 



After integrating out the k component, incorporating the contribution of Jfj° c , and 
using H20I21I52H . F(X j ,X j ) reduces to , 



F 3 +" (X j ,X j ) = 



li — v 



E 



exp 



J2 2eA 3 + eAA B (Y a ,6Y ) N?N? (36) 
1=1 /J 



2 

aJJEP 

= F a (X^Xj)-^ J2 B{Y a ,eY )MfM^ (37) 

a ,peP 

After integrating out the M-variables, and taking into account the observations fol- 
lowing {23 and (25, ^ a+n i x h X i) reduces to F 3 (Xj,Xj). This shows that in the 
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{Kj , Aj , Mj } representation the Aj are Gaussian random variables with positive mean 



<4> = 



2 {n-v) ^ 



(38) 



where 

and covariance 



<44> c = 



B[H k ,H{) 
1 



-r* 



(39) 



(40) 



// — v 

On the other hand, il38l4()[l together with ll-'ifill show that the n variables diffuse, while 
the a variables perform a ballistic motion, diffusing around straight trajectories. This 
means that for t 1 and with logarithmic precision the time dependent Lyapunov 
exponents coincide with the eigenvalues of 



2t 
1 



logT(a 2 (t)) 
-(exp^Kf))) 



— T 

t 



A (i') At' 



\ f At'Y,A k {t')r k 



(41) 



with T (exp (X j) — exp (r (X)) and r k = t (Hk). The operators all commute, so it is 
possible to write 14111 in block-matrix form. This means that the Lyapunov exponents 
are just the \ J dt' A k multiplied by the eigenvalues of the operators t^. 

Let us describe more precisely the statistics of the n variables. Expression ll37jl shows 
that the Mj are Gaussian random variables with zero mean and O (1) covariance. 
Expanding the Mj and Nj on the basis Y a one finds the explicit expression 



N a 

3 



exp 



-a 



-Ai 



Mf, 



(42) 



From the positiveness of a and of (Aj'), k = 1, . . . ,dimg, and passing to the con- 
tinuum limit, it follows that typical values of the N a (t) are exponentially small in 
t. Correspondingly, the n (t) variables become frozen after a time t — 0(1). The 
freezing of the n variables has been observed in Ref. ^H] and used to obtain analytic 
results about the statistics of the dissipation of a passive scalar field in Refs. |l()ll21j . 
Moreover, it provides a nice realization of Oseledec theorem about the stability of the 
Lyapunov basis in the t — > oo limit 32J. 
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7. Examples 



Let Eij denote a square matrix with entry 1 where the ith row and the jth column 
meet, all other entries being 0. Then 

EijEki = SjkEu, [Eij,Eki] = SjkEu — SuEkj 

The case of si (n, R). In this case gc is the algebra of n x n matrices with zero trace. 
Letting 

Hi = E u - E nn (1 < i < n - 1) , h c = CH i ( 43 ) 

i 

one has the direct decomposition [231 

gc = h c + CEij (44) 

i¥=3 

If H e he and e,; (H) (1 < i < n) are the diagonal elements of H, one has 

[H,E ij ] = (e i (H)-e j (H))E ij (45) 

The Killing form is 

B (X, Y) = 2nTr (XY) (46) 

The roots are 

e l - ej (1 <i,j < n) 

which are duals to the elements {Hi — Hj) £ hi?,. The compact real form for si (n, C) 
is u = su (n), the algebra of skew hermitian matrices; he are the diagonal matrices 
with null trace, a = hiR are the real matrices of he- Taking lexicographic ordering 
with respect to the basis a — e n (1 < i < n — 1) of the dual of Iir, i.e. assuming 
ej — e„ > ej — e n and consequently — ej > for i < j, one finds 

nc = CE^ = {strictly upper triangular matrices} 

i<j 

g = si (n, R) is the real form corresponding to the conjugation aZ = Z. The Cart an 
involution is OX = — X tr . One has k = u n g = so(n), p are the real symmetric 
matrices with null trace, a the diagonal matrices of p, n the real matrices of nc- The 
restricted roots coincide with the roots — ej. Letting Aj = X)fc=i AjHk and using 
il38l39l4r)l43l46j) . a simple computation gives: 

and 

Eq. (23J reproduces Eq. (39) in Ref. (23 letting ^ = l/4L>n, v = and Eq. (25) in 
Ref. j2U letting \x = l/2n 2 (n + 2) D, v = - (n + 1) /2n 2 (n + 2) D. 
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The case of sp (n, R). In this case gc is the algebra of complex matrices of the form 
X Y 



Z -X 1 



with Y = Y tr and Z = Z tr . Letting 



Hi = E u - E n+l . n+l (1 < i < n) , h c = J2 CH * ( 48 ) 

i=l 

one has the direct decomposition [23j 

gc = hp + ^ C (E n+ jj + E n+ j,j) + ^ C {Ej, n+ j + Ej, n+ j) 

i<j i<j 

+ J2c{E i , j + E n+j>n+i ) (49) 

Let be the linear form on h given by ej (Hi) — <5.y ; then 

[H, E n+itj + E n+jti ] = - (e. t (H) + ej (H)) (E n+i>j + E n+jji ) (i < j) 
[H,E itn+j + E jtn+i ] = (ei (H) + ej (H)) (E i;n+j + E jt n+i ) (i < j) 
[H, Eij — E n+ j. n+i ] = (e, t (H) — ej (H)) (Eij — E n+ j^ n+i ) (i ^ j) 

The Killing form is 

B (X, Y) = (2n + 2) Tr (XY) (50) 
The non zero roots are (± signs read independently) 

±2e t (1 < 1 < n) , ±e, ± e 3 (l<i<j< n) 

The roots 2ej are dual to the elements 4 ( ra 2 +1 ) Hi 6 hR. The compact real form is 



> = {(-? l):U + Ut=0,V = V*} 



u = sp (n 

Taking lexicographic ordering with respect to the basis ej (1 < i < n) of the dual of hR 
it is seen that the positive roots are + ej (i < j) and — ej (i < j). Consequently 

nc = ^ C (E i: n + j + Ejj l+i ) + 2^ ^ — E n+ j. n+i ) 

i<j i<j 

= | ( ^ ^ tr ^ , Y = Y tr , X strictly upper triangular 1 

g = sp (n, R) is the real form corresponding to the conjugation aZ = Z . One has 
k = sp (n) n so (2n) 

p = { ( J *l V * = * tr > * = Y = YtI > Y = Y 



a = {diagonal matrices of p} 
n = {real matrices of nc} 
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Restricted roots just coincide with roots. Letting Aj = J2k=i ^j^-k and using (I38l39l40l48lfi0|) , 
a simple computation gives 




) 



1 n-k+1 



k = 1, . . . , n 



2(ji-v) 2(»+l) ' 



and 




1 



6 kl 



4 (n + 1) (/i - i/) 
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A. Lie-Cartan theory 

I recall here some results from the theory of Lie algebras 23, 24, 27, 28 for convenience 
and in order to establish the notations, which partially parallel those of Ref. [231 • 

Any complex Lie algebra gc can be seen as a real Lie algebra endowed with a 
complex structure. A real form of a complex Lie algebra gc is a real Lie algebra g 
such that gc = g + zg (direct sum of vector spaces) where g and ig are identified with 
subalgebras of gc- Thus gc is isomorphic to the complexification of g. (E.g., si (n, R) 
is a real form of si (n, C)). There is a 1-1 correspondence between the real forms of a 
Lie algebra and the conjugations a : X + iY — > X — iY (X, Y G g). (E.g., si (n, R) 
corresponds to the conjugation a : Z — > Z of si (n, C)). A compact Lie algebra is a 
real Lie algebra which is isomorphic to the Lie algebra of a compact Lie group. Every 
semisimple Lie algebra gc over C has a compact real form u, which is unique (up to 
automorphisms; e.g., su (n) is the compact real form of si (n, C); it corresponds to the 
conjugation t : Z — > —Z'). 

A Cartan subalgebra he of a complex Lie algebra gc is a maximal Abelian subal- 
gebra such that the endomorphism adiJ : X — > [H,X] is semisimple (i.e., gc can be 
written as a direct sum of irreducible invariant subspaces of adiJ for all H £ he)- 
Every semisimple Lie algebra over C has a Cartan subalgebra which is unique (up to 
automorphisms; e.g., diagonal matrices with null trace form a Cartan subalgebra of 



On each real and complex Lie algebra one can define the Killing form B (X, Y) = 
Tr (ad X ad Y), where Tr denotes the trace of the vector space endomorphism. (E.g., 
B (X, Y) = 2nTr XY for X, Y 6 si (n, C)). Every semisimple Lie algebra over C has 
a real form u such that B is strictly negative definite on u (it then follows that u is 
compact, cf. above; e.g. TrX 2 = -TrXX 1 < for X 6 su (n), X ^ 0). 



si (n,C)). 
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A complex Lie algebra gc over C can be decomposed in the direct sum of the 
simultaneous eigenspaces (root subspaces) 

gc = {X £ gc : ad(JT) X = a (H) X for all H £ h c } 

where the linear function a : he — > C is called a root (of gc with respect to he) if 
g c 7^ {0}. Let A denote the set of nonzero roots. One has the root decomposition 

gc = h c + g c (direct sum) (51) 

with the properties: i) dimg c = 1 for each a £ A; ii) g c and g c are orthogonal under 
B if a, (3 are any two roots and (3 ^ —a; Hi) B is non degenerate on he and for each 
linear function a on he there exists a unique H a £ be such that 

a(H) = B(H,H a ) (52) 

for all H £ h c ; iv) ii a £ A, then -a £ A and [g^gc"] = CH a and a(H a ) ± 0. 
Let hR = X)aeA R-ffct- Then by m,), A can be seen as a subset of the dual space h R . 
Given any basis on Iir one can introduce a lexicographic order on h R and thus on A. 
This gives the concept of positive and negative roots. A positive root is called simple 
if it cannot be written as the sum of two positive roots. There are exactly dim ha 
simple roots aj and each a £ A can be written as J2 n j a j with nj integers which are 
either all positive or all negative. 

For each a £ A one can choose X a £ g c such that for all a, j3 £ A: 

[X a ,X- a \ = H a 

[H, X a ] = a (H) X a for H £ h c 

[X a ,Xfs] = if a ^ -[3 and a + & A 

[X a , X p ] = N a ^X a+fj if a + [3 £ A, with N aJj = -N_ ay _ 

B(X a ,Xp) = 8 a +f3fl 

On the real algebra u generated by the elements iH a , X a — X- a , i (X a + X- a ), a £ A, 
the form B is strictly negative definite; so, u is a compact real form of gc (c/. above). 
If r is the conjugation of gc with respect to u, rX a = —X- a for a £ A. The set {X a } 
is called a Weyl basis. 

A Cartan decomposition of a semisimple Lie algebra g over R is a direct sum g = k+p 
in a subalgebra k and a vector subspace p such that the mapping 9: K + P^K — P 
(K £ k, P £ p) is an automorphism of g and the form —B(X, 9Y) is strictly positive 
definite. 6 is called a Cartan involution. Each semisimple Lie algebra over R has 
a Cartan decomposition. If gc is the complexification of g, a is the conjugation of 
gc with respect to g, one can always choose a real form u such that cm C u, so 
that k = uHg, p = iu n g and u = k + ip; u is then a compact real form of gc- 
(E.g., for g = si (n, R) the Cartan decomposition is the decomposition of a matrix 



19 



into its antisymmetric and symmetric part; in this case OX = — tT X, k = so(n), 
p = {symmetric matrices} and k + ip = su (n)). If gc is a semisimple Lie algebra 
over C, u its compact real form, then gc = u + iu is a Cartan decomposition of gc- 
(E.g., for gc = si (n, C) it is the decomposition of a matrix into its skew-hermitian 
and hermitian part). 
Since 

6 = 1 on k and = -lonp (53) 

then, for all 1,7 eg, 

0&d(X)Y = a,d(0X)0Y (54) 
6>Ad(exppf))y = Ad (exp (OX)) BY (55) 

where Ad is the adjoint representation of G on g (e.g. Ad(g) X — gXg^ 1 if G, g are 
linear), so that Ad (exp (Xj) = exp (ad (X)). 

Let h be a maximal Abelian subalgebra of g, with a = h n p as large as possible, 
and let ai = h n k. Then Pir = a + ia% and he = hR + iliR is a Cartan subalgebra of 
gc- (E.g., for g = si (n, R), ai = {0} and a = hR are the real diagonal matrices with 
null trace) . Let P be the set of positive roots that do not vanish identically on a (for 
si (n, R) this just coincides with the set of positive roots). Let 

nc = ^Sc> n = n c ng 

Then nc and n are nilpotent Lie algebras and one has the direct vector space sum 

g = k + a + n (56) 

(Iwasawa decomposition). From ii) it follows that a is orthogonal with respect to B to 
both k and n, and that n is orthogonal to itself. The decomposition 15 fill exponentiates 
to 

G = KAN (57) 
where K = exp k, A = exp a, N = exp n, i.e., the mapping 

(k, a, n) — > kan (58) 

is a diffeomorphism of K x Ax N onto G. (E.g., for si (n, TV), k is the rotation algebra, a 
is the algebra of diagonal matrices with null trace, n is the algebra of strictly triangular 
matrices; at the group level this decomposition can be naturally obtained through the 
Gram-Schmidt orthogonalization process). 

To the root decomposition l|5ll of gc there corresponds a restricted root decompo- 
sition 

g = a + m+ g"' 

q'GA' 
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where m are the elements of k which commute with all the elements of a, A' is the set 
of nonzero restrictions a' of a to a and 

g Q ' = {X G g : ad(if) X = a' (H)Xfov all H G a} 

Since different roots a can have the same restriction a' to a, the multiplicity 

m a i — dimg" (59) 

can be higher than 1 for some algebras. The Cartan involution changes the sign of the 
restricted roots: 

0g a ' = g- a ' (60) 
The algebra k can therefore be decomposed as 

k=m + E (s a '+H a ') (61) 

A Weyl basis {X a } for gc can be projected on g by 

Y a = l - (X a + aX a ) (62) 

so that g a =J2 R-^a where the sum runs over all a £ A whose restriction to a is a'. 
On a Lie group one can define the (right invariant) Haar measure 

dimG 



d.9 = f\ 9 1 Sig 



i=i 

where Sig = g SiX and the SiX are a basis of 1-forms on g. Semisimple Lie groups are 
unimodular, i.e. such that |det Ad (g)| = 1 for each g € G, so that right invariant mea- 
sures are also left invariant. The Iwasawa decomposition induces the decomposition 

dg = j (k, a, n) dk da dn (63) 

where j is a Jacobian factor. The groups G, K, A, N are all unimodular, therefore 
both sides of ll63)l are invariant under the right and left translations g — * k^g, g — > gn^ 
with &r G K and til G N, so that j (k, a, n) = j (a). On the other hand, the l.h.s of 
il63H is also invariant under g — > ga^ 1 , which corresponds to kan — > fee - (ana -1 ) (with 
e the identity element), so that 

. , , . , , d (ana' 1 ) 

J ( a ) = J (e) 



d(n) 

= j(e)det(Ad(o)| n ) 

= j (e) exp [Tr ad (exp -1 a) | ] 



j (e) exp 



a (exp 1 a) 
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If the measures on G, K, A, N are conveniently normalized we can assume that 
j (e) = 1, so 

dg — exp [2p (exp" 1 a)] dfc da dri (64) 

with 

P = 2 H 01 ( 65 ) 

Let r be a skew-hermitian representation of the compact form u on a finite dimen- 
sional Hilbert space V . By analytically extending r to the complexification gc = u + iu 
and then restricting to the real form g = k + p one obtains a representation of the 
real algebra g (Weyl unitary trick). If K 6 k, P G p, then r (K) and r (iP) are 
skew-hermitian, while r (P) is hermitian. In other words, 

r f (X) = -t{6X) 

Denoting T (exp (JQ) = exp (t (X)), this property exponentiates (23 p. 144] to 

Tt (g) = T- 1 (Qg) (66) 
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